Quasi-adiabatic Continuation of Quantum States: The Stability of Topological 
Ground State Degeneracy and Emergent Gauge Invariance 
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We define for quantum many-body systems a quasi-adiabatic continuation of quantum states. 
The continuation is valid when the Hamiltonian has a gap, or else has a sufficiently small low- 
energy density of states, and thus is away from a quantum phase transition. This continuation takes 
local operators into local operators, while approximately preserving the ground state expectation 
values. We apply this continuation to the problem of gauge theories coupled to matter, and propose 
a new distinction, perimeter law versus "zero law" to identify confinement. We also apply the 
continuation to local bosonic models with emergent gauge theories. We show that local gauge 
invariance is topological and cannot be broken by any local perturbations in the bosonic models in 
either continuous or discrete gauge groups. We show that the ground state degeneracy in emergent 
discrete gauge theories is a robust property of the bosonic model, and we argue that the robustness 
of local gauge invariance in the continuous case protects the gapless gauge boson. 



INTRODUCTION 

Traditionally, gauge theory was described as a theory 
of a vector field that has a local gauge symmetry. In a 
Lagrangian framework, this means Lia^ + d^cf)) = L(a^). 
This gauge symmetry was believed to be the defining 
property of a gauge theory. It was believed that the 
gauge symmetry protects the gapless gauge boson for 
continuous gauge groups and the topological ground state 
degeneracy on compact space for discrete gauge groups. 
Even the slightest gauge symmetry breaking, such as a 
gauge potential term (a M ) 2 in the Lagrangian, gives a 
finite mass to the gauge boson or lifts the topological de- 
generacy. In this case we can no longer regard the theory 
as a gauge theory at low energies. 

However, the above standard picture for gauge sym- 
metry protecting gapless gauge bosons and ground state 
degeneracy is very formal, since the gauge symmetry is 
not really a symmetry. Within the Hamiltonian formula- 
tion of gauge theories, the gauge transformation is simply 
a transformation between different labels that label the 
same physical state. It is a do-nothing transformation. 
It is very different from the usual symmetry transforma- 
tion that transforms a physical state to a different phys- 
ical state. Therefore, it is not clear what is the essence 
of gauge theory and gauge symmetry. 

In last 15 years, it was shown, with increasing rigor, 
that deconfined gauge theories can emerge from certain 
local bosonic lattice models. [1-16] If the emergent gauge 
theory has a continuous gauge group, the local bosonic 
model will have gapless excitations that behave just like 
gauge bosons. [8, 11, 14-16] If the emergent gauge theory 
has a discrete gauge group or has a Chern- Simons term, 
the local bosonic model will have degenerate ground 
states on compact space. [4, 17, 18] This raises a physi- 
cal question: what protects the gapless gauge bosons and 



the ground state degeneracy? According to the standard 
picture for gauge theory, those properties are protected 
by gauge symmetry, but from the point of view of local 
bosonic model, what is this "gauge symmetry"? How 
does "gauge symmetry" emerge at low energies? 

A close examination of those local bosonic models 
with emergent gauge theory reveals that the emer- 
gence of gauge theory is intimately connected to string 
condensation. [11, 13] The "gauge symmetry" is related to 
the integrity of the strings. If the strings are unbreakable, 
one can show that the low energy states are gauge invari- 
ant. However, in general, the strings in the boson model 
are not perfectly well defined. Strings may break up mo- 
mentarily and rejoin. One may wonder if this means 
that the "gauge symmetry" becomes approximate. We 
know that a theory that loses its "gauge symmetry" even 
slightly no longer behaves like a gauge theory at low en- 
ergies. This seems to suggest that breakable strings will 
give gauge bosons a mass gap or lift the ground state 
degeneracy. 

On the other hand, it was believed that the gaplessness 
of the gauge bosons and the degeneracy of the ground 
states in those bosonic models are topological and are ro- 
bust against any local perturbations. [4, 17, 18] A formal 
argument goes as the follows (see, for example, Ref. [19], 
page 393 and page 435). We first derive the low energy 
effective gauge theory of the bosonic model. We then ar- 
gue that any generic perturbation of the bosonic model 
cannot generate any terms that break the "gauge symme- 
try" in the low energy effective gauge theory. Therefore, 
all the properties protected by the "gauge symmetry" 
are robust against arbitrary perturbations of the bosonic 
model. 

We see that to understand why "gauge symmetry" re- 
mains exact even for a generic boson model with virtually 
breakable strings is vital in our understanding why the 



degeneracy of the ground states is protected even when 
the original boson model has no symmetry, and why the 
gaplessness of the gauge bosons is protected even when 
the original boson model has only translation symmetry. 
We would like to address some of these issues in this pa- 
per. 

In the next section, we start by introducing some 
bosonic models with emergent gauge theories, and raise 
the question of the ability of topological order and gauge 
symmetry to survive local perturbations of the Hamilto- 
nian, as well as proposing the zero law to identify de- 
confined gauge theories. We then introduce in the fol- 
lowing section a quasi-adiabatic continuation which en- 
ables us to identify appropriately dressed operators for 
the perturbed Hamiltonian, such that the dressed oper- 
ator has almost the same ground state expectation value 
for the perturbed Hamiltonian as the original operator 
had for the original Hamiltonian. The final section il- 
lustrates the application of the continuation via a series 
of examples, beginning with local models without emer- 
gent gauge structure, such as quantum Ising models, and 
then going on to emergent gauge theories. Many of these 
systems are theories for which the unperturbed Hamil- 
tonian has a ground state degeneracy and then a gap to 
the rest of the spectrum, with no local operators con- 
necting the degenerate ground states. We are then able 
to show that, so long as the gap to the rest of the spec- 
trum remains open, the splitting between the low energy 
states of the perturbed Hamiltonian remains exponen- 
tially small. This is illustrated in the case of the ferro- 
magnetic quantum Ising model; in the case of the frac- 
tional quantum Hall effect where we are able to extend 
results on the insensitivity of the topological degeneracy 
to disorder [17, 18]; and in the case of an emergent 
gauge theory. The case of gapless theories is also dis- 
cussed; it will turn out that gauge symmetry is much 
more robust in compact than in non-compact theories. 



SIMPLE LOCAL BOSONIC MODELS WITH 
EMERGENT GAUGE THEORIES 

To make the above discussion more concrete, in this 
section we are going to discuss two simple bosonic models 
that have emergent Z 2 and U(l) gauge theories respec- 
tively. 



A bosonic model with emergent Z2 gauge theory 

The first bosonic model is a spin- 1/2 model on a d- 
dimensional cubic lattice. [20-22] The spins live on the 
links labeled by i. The Hamiltonian is given by 



H = uY,(l-W I )-gJ2\ II 



(1) 




FIG. 1: A closed-string state. The up-spins are represented 
by open dots and the down-spin by filled dots. 
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FIG. 2: The energy levels of N spin-1/2 spins described by 
H . 



where I labels the vertices and p the squares of the lat- 
tice, and Wj is given by 



Wi 



n * 

legs of I 



(2) 



The legs of a vertex are the links that connect to the 
vertex and the edges of a square are the four links around 
the square. a x ' v,z are the Pauli matrices. 

When U is very large, the low energy sector of the 
model is formed by closed string states. What are the 
closed string states? First, the state with all spins up 
is defined as the no string state. A closed string state 
is a state where the down spins form closed loops (see 
Fig. 1). We note that the close-string states are the 
states that satisfy Wj | closed strings) = | closed strings), 
Since [Wj, H ] = 0, the closed-string states and the open- 
string states [33] do not mix. This allows us to plot the 
spectrum of Ho separately for closed-string states and 
open-string states in Fig. 2. 

When g > 0, the ground state |$q) of the model is 
the equal weight superposition of all closed string states. 
Such a state is called a closed-string condensed state since 
the closed-string creation operator 5(C c ioscd) has a non- 
zero expectation value [23] 



($ |S(Ccloscd)|$0) = 1 



(3) 



regardless of the size of the string C c i OS ed- Here a string 
creation operator that creates a string C is given by 
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FIG. 3: A dual closed-string. 



We can show (3) to be true by noting that 

[S^Cclosed ),H ]=0. (4) 

What is the physical character of the closed-string con- 
densed state? It is believed that the closed-string con- 
densed state contain a new kind of order - topological or- 
der - which cannot be characterized by symmetry break- 
ing and long range order. [13, 17] So we need a new way 
to characterize such an order. One way to characterize 
the topological order is through the robust ground state 
degeneracy on torus. [4, 18] The closed-string condensed 
state in our spin-1/2 model is characterized by a four-fold 
ground state degeneracy. 

To show the four-fold ground state degeneracy, we 
would like to first point out that the S(C c \ ose ) are not 
the only closed string operators that commute with the 
Hamiltonian. We can define dual string operators that 
also commute with the Hamiltonian. [22] 

While strings are formed by segments that connect 
nearest neighbor vertices, dual strings are formed by 
segments that connect the centers of nearest neighbor 
squares (see Fig. 3). A dual string operator S(C) for a 
dual string C is defined as 

S(C) = 1] of. (5) 

i cross C 

One can check that 

[S(C c i osed ), H ] = 0, (6) 
which implies that the dual closed-strings also condense 

($0|S(4lo S ed)|$0) = 1- (7) 

Now we are ready to show that the ground states 
of H on torus have at least four-fold degeneracy. Let 
C x and C y (C x and C y ) be the closed (dual) strings 
that wrap around the torus once in x and y direc- 
tions. We find that the four large-string operators 
(S(C X ), S(C y ), S(C X ), S(C y j) commute with each other 
except 

{S(C x ),S(C v )} = 0, {S(C x ),S(C v )} = 0. (8) 



The above algebra has only one four-dimensional irre- 
ducible representation. Since the large-closed-string op- 
erators (S(C x ),S(Cy),S(C x ),S(C y )) act within the de- 
generate ground states, the ground state degeneracy 
must be a multiple of four. 

When U = oo, the model becomes the gauge 
theory. [20] One way to see this is to note that the only 
states with finite energies are closed-string states and 
closed-string states are gauge invariant states under local 
Zi gauge transformations. The local Z^ gauge transfor- 
mations are generated by the unitary operators Wi in 
(2): 

Wi | closed strings) = | closed strings) for any I. (9) 
A generic Zi gauge transformation is given by 

J|(Wj) ni | closed strings) = |closed strings). (10) 

In the gauge theory language the closed-string operator 
'S'(Ccioscd) turns out to be the Wilson-loop operator, [20, 
24] which is gauge invariant 

WiS{C c \ osc d)W] = S(C C 

loscd ) • 

(3) implies that the expectation values of the Wilson 
loop satisfies the perimeter law 

($o|S(Cdoscd)|$o) ~ e~^ c ^ (11) 

with zero coefficient a = 0. Here |C c i OS cd| is the length 
of the string C c i osc . In this case we will call (3) the zero 
law. 

The perimeter law indicates that the Z2 gauge theory 
is in the deconfined phase. The Z2 deconfined phase has 
four nearly degenerate ground states on the torus which is 
consistent with our previous direct calculation on the spin 
model. The energy separation between the four nearly 
degenerate ground states is of order e~ L ^ where L is the 
linear size of the torus and £ a finite length scale. 

Now let us add a term 

i 

to our spin model Hamiltonian Hq and assume U is finite. 
S(C X ) and S(C y ) no longer commute with the modified 
Hamiltonian Hq + Hi. So it is not clear if Hq + Hi still 
has four degenerate ground states on a torus. 

To understand the properties of modified spin system 
Hq + Hi, we note that Hq + Hi still does not mix the 
closed-string and open-string states. So if U ^> g, Ji, the 
low lying states are still closed string states. One can 
check that if we restrict Hq + Hi to the close-string sub- 
space, the system is identical to a pure lattice Z2 gauge 
theory. A pure Z2 gauge theory has two phases in d+1 
dimensions if d > 1: a deconfined phase where the ex- 
pectation value of the Wilson loop satisfies the perimeter 
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FIG. 4: The energy levels of N spin-1/2 spins described by 
H a + H! + H 2 , assuming |Ji|, |J 2 | < |g| < U < JV|g|. 



law, and a confined phase where the expectation value of 
the Wilson loop satisfies the area law: 



($o|S(C closcd )|<I>o}~e- 



-lA{C c 



i) 



(12) 



where A(C c iosod) is the area enclosed by the loop C c iosod 
and 7 > 0. So we expect our spin model H + Hi also 
has two phases. If | Ji | <C \g\, the ground state of the spin 
model is filled with large closed strings which correspond 
to the Zi deconfined phase. If J\ 3> |.g|, the ground state 
of the spin model hardly has any strings (i.e. almost all 
spins point up) which correspond to deconfined phase. 

Base on this picture, we guess that when |Ji| <C 
the ground states of H + Hi are (nearly) four-fold de- 
generate, while when |Ji| » \g\ the ground state is not 
degenerate. But the result for | Ji| <C \g\ clearly is just a 
guess. Can we provide a more rigorous proof? 

The situation gets even more complicated if we add 
another term 



H, 



Such a term can mix the closed-string and open-string 
states. The low energy sector of H + H i + H 2 is not 
formed by simple closed-string states, although the mix- 
ing with the open string states may be small for small J 2 
(see Fig. 4). It appears that Hi breaks the Z 2 "gauge 
symmetry" since the low energy states no longer satisfy 
(9). 

Also when J 2 ^ and U < oo, the expectation value 
of a closed string operators satisfies the perimeter law in 
both \Ji\ <C \g\ limit and J\ 3> \g\ limit, which give no 
sign of two phases. All of those suggest that the four-fold 
ground state degeneracy is lifted by a finite J 2 . 

But this suggestion is incorrect. Ref. [4] argues that 
any local perturbations of H cannot break the Z 2 "gauge 
symmetry" . As a result, Ho + Hi + H 2 will have four 
nearly degenerate ground states as long as Ji and J 2 are 
not too large. Such a phase contains a non-trivial topo- 
logical order. However, Ref. [4] only provides a formal 
argument. A more rigorous understanding is needed. 
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FIG. 5: A likely quantum phase diagram for the spin-1/2 
system Ho + Hi + H2 ■ The deconfined phase is characterized 
by four nearly degenerate ground states on torus with energy 
splitting of order e~ L ^ where L is the linear size of the torus 
and £ a length scale. The confined phase is characterized 
by a non degenerate ground state on torus. In general, the 
confined phase and the deconfined phase are distinguished by 
the zero law and the perimeter law of certain loop operators. 



Certainly, large Ji and J 2 will polarize the spins and 
lift the ground state degeneracy. Such a phase has a 
trivial topological order. This suggests a phase diagram 
in Fig. 5. We note that the two phases in Fig. 5 have the 
same symmetry and are distinguished only by topological 
orders. 

From the phase diagram, we see that when J 2 7^ 
the perimeter/area laws of the closed string operators 
(or Wilson loop) cannot determine if the ground state is 
topologically ordered or not. Thus for a generic bosonic 
model, perimeter/area laws of the closed string opera- 
tors is not the proper way to test if the ground state 
has closed-string condensation (or non-trivial topological 
order). 

In the examples section, we will show that the topolog- 
ical phase is characterized by dressed closed-string oper- 
ators which satisfy the zero law. The trivial phase does 
not contain any such closed string operator. The closed 
string operators in the trivial phase all satisfy the perime- 
ter law. So it is the zcro/pcrimetcr laws that distinguish 
topological/trivial (deconfined/confined) phases, instead 
of perimeter /area laws. 

For J 1 = J r 2 = 0, the system has an exact four- fold 
degeneracy of the ground state on the torus, followed by 
a gap of order g to the next lowest state. The closed- 
string operators satisfy the zero law. In the examples 
section, we will show that, for small but non-zero Ji and 
J 2 , a deformed or dressed closed-string operator can still 
satisfy the zero law. The zero law of the dressed string 
operator allows us to show the four-fold degeneracy of the 
ground states in the small Ji, 2 limit. More precisely, we 
assume that the gap, from the four lowest states to the 
rest of the spectrum, remains open, and then we show the 
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FIG. 6: The empty dots represent rotors with L z — - the no 
string states. A closed string is formed by a loop connecting 
neighboring vertices. A closed-string state is obtained by al- 
ternately increase or decrease L z by 1 along the closed string. 
The filled dots represent rotors with L z = ±1. The arrows on 
the links all point from even vertices to odd vertices. 



four-fold ground state degeneracy up to an exponentially 
small splitting. We will also show that, in the small Jx >2 
limit, the low energy sector of the model is still formed by 
Z 2 gauge invariant states. However, for J 2 ^ 0, the gauge 
invariance is under a deformed Z 2 gauge transformation. 
So in this sense, none of the small perturbations in the 
spin model can break the "gauge symmetry" in the low 
energy effective gauge theory. 

A bosonic model with emergent (7(1) gauge theory 

In our second bosonic model, we consider rotors on 
the links of d-dimensional cubic lattice. A rotor can be 
viewed as a particle moving on a circle. The position of 
the particle is given by an angle 9, and the angular mo- 
mentum of the particle by L z — —idg. The Hamiltonian 
of the rotor model is given by 

ffrotor = U Q 2 I - 9 J2 ( B P + ^- C -) 
I P 

+ J 1 5>f) 2 + J 2 ]>>++L7) (13) 



Bp — L^L 2 L^L 4 , 



(-) 



legs of I 



where I labels the vertices, i labels the links and p labels 
the squares of the cubic lattice. 1, 2, 3, 4 label the four 
links that form the edges of the square p. L + = e 10 is 
the raising operator of L z , L~ — (L + y , and (— ) I = 1 for 
the even vertices and. (— ) 7 = —1 for the odd vertices. 

When g = J\ = J 2 = and U > 0, the ground states 
are highly degenerate and form a low energy subspace. 
One of the ground states is the state with L| = for ev- 
ery rotor. Other ground states can be constructed from 
the first ground state by drawing an loop in the cubic lat- 
tice, and then alternately increasing or decreasing L z by 
1 along the loop (see Fig. 6). The sum J2 lcgs of j L\ van- 
ishes on every vertices for such a closed-string state. Such 



a process can be repeated to construct all of the degen- 
erate ground states. We see that the degenerate ground 
states are formed by loops, or more precisely string-nets, 
since loops can overlap and form branched strings. 

The Ji-term gives the strings in the degenerate ground 
states a fine energy and represents string tension. The B p 
operator creates a small loop of closed string or deforms 
the existing strings. Thus the g-term generates string 
"hopping" or string fluctuations. 

The Ji^-term and the g-term lift the degeneracy of 
the ground states. In the U ^> J\ S> max(|g|, | J 2 |) limit, 
the true ground state correspond to a state with almost 
no strings (i.e., a state with L\ — on every link). The 
excitations above such a state have finite energy gaps. In 
the U ^> \g\ S> Ji,2 > limit, the true ground state is a 
superposition of many large closed strings [11, 19]. Such 
a state is a string-net condensed state. 

When J 2 = 0, the Hamiltonian H ro t OI does not mix 
the closed-string states and open-string states. When 
restricted to the closed-string subspace, H rotOT is iden- 
tical to the Hamiltonian of lattice U(l) gauge theory. 
The closed-string states can be viewed as gauge invari- 
ant states since they are invariant under local U(l) gauge 
transformations 

e i<pQi | closed-string) = | closed-string), for any J. (14) 

A general U(l) gauge transformation is generated by 
e iY.i<f>iQi _ 

In the limit \J\\ ^> \g\, the lattice f/(l) gauge theory 
is in the strong coupling limit and is in a confined phase. 
In the limit \J\\ -C |<?|, the lattice U{\) gauge theory is 
in the weak coupling limit and has gapless U(l) gauge 
bosons as its only low lying excitations if d > 2. So when 
J 2 = and when U 3> \g\ ^> | Ji|, the rotor model H Iotor 
contains emergent gapless U(l) gauge bosons. [11, 14-16, 
19] 

When J 2 7^ 0, H rotor mixes the closed-string and open- 
string states. The low energy states are no longer pure 
closed-string states and is not invariant under the local 
U(l) gauge transformation (14). It appears that a non- 
zero J 2 will break the (7(1) gauge symmetry. We may 
conclude that even a small J 2 will give the [7(1) gauge 
boson a gap and the rotor model H I0t0 r cease to have 
emergent [7(1) gauge bosons at low energies. In the sec- 
tion of examples, we will show that this line of argument 
is incorrect. For a small J 2 (or any other small perturba- 
tions to i/rotor), we can define deformed local [7(1) gauge 
transformations so that the low energy states of H rotm are 
invariant under the deformed local [7(1) gauge transfor- 
mations. Thus the local "[7(1) gauge symmetry" cannot 
be broken by any small perturbations if U ^> |^i,2|> 
Thus far, we will prove these results for [7(1) theories; 
we do not prove, but strongly conjecture, that the stabil- 
ity of the [7(1) gauge symmetry protects the gaplessness 
of the gauge boson. As a result, no translation invari- 
ant perturbation can give the gapless [7(1) gauge boson 
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a mass gap. The gaplessness of the emergent U(l) gauge 
boson is topologically stable! 



QUASI- ADIABATIC CONTINUATION 

In this section we define the quasi-adiabatic continua- 
tion. We consider a family of Hamiltonians, H s , depend- 
ing on a continuous parameter s, where we wish to define 
a continuation from s = to s = 1. We denote eigen- 
states of H s by ^ a ,s, where a state ^ a ,o evolves into state 
ty a ,s under an adiabatic evolution of H s - In the event of 
a level crossing as a function of s, any arbitrary continu- 
ation of the states through the level crossing is allowed. 

Let us begin with some motivation and discussion. 
We define the unitary operator V(s) = J2 a \^a,s)(^a,o\- 
Then, for any operator O, we could define a correspond- 
ing dressed operator by O a diab(s) = V(s)OV(s)^ so that 
Oadiab(s) would have exactly the same expectation value 
in state \I/ ajS as O does in state "J/a.o- Using such a defi- 
nition of dressed operators, we can show that the dressed 
string operators and the dressed gauge transformations 
will have the same properties in the deformed model as 
the bare string operators and the bare gauge transfor- 
mations in the exactly soluble model. However, this def- 
inition of a dressed operator would not suffice for our 
purposes at all! In particular, we do not have any reason 
to believe that the resulting O a diab(s) would still be a 
local operator. 

Indeed, even if O only involves operators on a few sites, 
in general, the continued O a diab(s) will contain operators 
on every site of the system. However, in this section, we 
will show that, under certain conditions, we can continue 
O into an local operator that acts only on a finite number 
of sites. 

To state the result more precisely, let us first assume 
that H s has a gap AE separating a low energy sector and 
a high energy sector for all < s < 1. For any operator 
O which acts on a set of sites So, we will construct an 
approximate dressed operator O(s) which only acts on 
sites within a distance I from the sites So- Then, there 
is a unitary matrix Qo(s) which acts only within the low 
energy sector of H s , such that 

(*iow, s |Qo(s) t O(s)g (s)-O adiab ( S )|* low , s ) < N Soi e~ 1 ^ 

(15) 

in the large I limit for a certain fixed length scale £, where 
Ns o l is the number of sites within distance I of a site in 
So- This number grows only as a power of Z, so is easily 
overwhelmed by the exponential. The state l^iow.s) is 
any state in the low energy sector of H s ■ 

If H s has no gap between the low energy and high 
energy sectors, we will define a density of states p(E) in 
the high energy sector. If p(E) is bounded by p(E) < 



E a 1 , then the bound in (15) is weakened to 

(*iow, s |Qo(s) t O( S )Qo(s)-O adiab (s)|* 1 ow,s) < Nso.l 1 -^ 2 - 

(16) 

For a local operator O, Ns Q t grows as l d , where d is the 
dimension of the lattice, and soifrf+1 — a/2 < then 
the error decays for large I. For a string-like operator, 
Ns ,i i s proportional to Z d_1 |0|, where \0\ is the length 
of the strength, for I < \0\, and is proportional to l d for 
I > \0\. (A more rigorous statement of the results will 
be given later.) 

The key in obtaining the above result is to adopt a 
different definition of the dressed operator by O(s) = 
V(s)OV(sY , where the unitary operator V(s) is defined 
following Eq. (17). Physically, the definition of V corre- 
sponds to adiabatically changing the Hamiltonian from 
Ho to H s , while the definition (17) of V(s) corresponds 
to a quasi-adiabatic change of the Hamiltonian. So V(s) 
can be viewed as an approximation of V(s). The opera- 
tor V(s) will be chosen to achieve the goal of defining 
dressed operators O(s) which have approximately the 
same ground state expectation values in the perturbed 
Hamiltonian H s as the original operators did in the un- 
perturbed Hamiltonian, while preserving the locality of 
the operators. 

The locality of the operators is very important. Sup- 
pose, for example that the ground state of Ho has long- 
range correlations. That is, there exist two local oper- 
ators, Oi , O2 which are correlated even though the two 
operators are far separated from each other in space. For 
example, if this is a spin system with long-range spin 
correlations these may be spin operators acting on two 
different sites which are far separated from each other. 
Then, however, the operators 0\(s),02(s) will also be 
correlated in the ground state of H s , and since the oper- 
ators remain local under the continuation this implies the 
existence of long-range correlations in the ground state 
of Hs- 

The major result will be an explicit definition of V(s) 
in terms of derivatives of the Hamiltonian which accom- 
plishes these goals. In order to show that ground state 
expectation values remain approximately unchanged un- 
der the continuation, we will make some assumptions on 
the existence of a gap, though some extensions to gap- 
less systems with sufficiently small low-energy density of 
states will be discussed. The dressed operators remain 
local under this evolution: local operators are spread out 
over a length scale of order the correlation length in a 
gapped system, while string-like operators are spread out 
over a length scale which grows only logarithmically with 
the string length. Detailed proofs will be given in an Ap- 
pendix. 

The importance of the low energy density of states 
in what follows can be understood physically by anal- 
ogy to another continuation that should be much more 
familiar, namely Fermi liquid theory. As discussed by 
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Anderson[25], the correct way to think of Fermi liquid 
theory is to think of starting with a non-interacting sys- 
tem and turning the interactions on slowly but not in- 
finitely slowly; that is, physically exactly the same proce- 
dure we imagine here. Anderson's discussion of how fast 
the interactions need to be turned on is based on consid- 
erations of the quasi-particle states, and corresponds very 
closely to our two criteria: maintaining both the expec- 
tation values and the locality of the operators. To show 
that this is possible, the analysis in the Fermi liquid case 
relies on the low density of particle-hole excitations near 
the Fermi surface. Here, we rely on something similar, 
namely a low density of states at low energy. Our contin- 
uation is very general, and thus valid for a much wider 
range of systems than Fermi systems, but this generality 
can in some cases limit what we can prove on specific 
systems. 



Definition of System 

We consider a family of Hamiltonians, TL S1 which obey 
the finite-range conditions [26, 27]: TL S = J^i^li where 
letters i, j, ... labels different lattice sites of the system; 
each ri\ acts only on sites j with d(i,j) < R where R 
is the interaction range and d(i,j) is some metric on the 
lattice; and | \H\ \ \ < J for some constant J for all i, s. We 
further assume that 1 1 d s Hl \ \ < K for some constant K for 
all i, s. It is possible to slightly weaken the finite-range 
conditions and consider exponentially decaying interac- 
tions as well [26]. 

To define the concept of "density of states", let us fo- 
cus on a state ^low.s that we continue. We will consider 
a 'J'low^ which is a ground state of Ti Sl but in general we 
could continue any eigenstate of fi s . The state ^ a , s has 
energy E a > compared to ^low.s if a is not in the low 
energy sector. We assume that, for < s < 1, the den- 
sity of states of the local operator u\ = d s W s is bounded 
as follows. We assume that there is a D{E) such that 
Y.aeM gh ,\E a \<E\^aM\^., s )\ 2 < D(E)\\u% where 
||...|| denotes the operator norm. The sum is restricted to 
states <J/ a in the high energy sector. Note that D(E) < 1 
for all E. We define 8eD{E) = p(E) to be the density 
of states at energy E produced by operator u l s . 

There are a number of systems for which this density 
of states bound is relevant. For a discrete gauge theory, 
there will be a set of topological excitations below a gap. 
These topological excitations form the low energy sector, 
and thus D(E) vanishes below the gap. For a transverse- 



<i,j> °i a j 



Ba- 



the 



field Ising system, W = JJ2 

ferromagnetic phase, there will again be two almost de- 
generate states below a finite gap. At B = 0, these states 
correspond to symmetric and anti-symmetric combina- 
tions of all spins up or all spins down. 

The bound D(E) implies a locality bound[27]. Sup- 
pose some operator O acts only on some set of sites, 



So- We define the distance between a site j and the 
operator O to be equal to d{j,0) = mim e s (d(i, j)), 
with the minimum is taken over sites i e So, and we 
define the distance between two operators Oi , O2 to be 
d(0 1 ,0 2 ) = mm ieSoi je s 02 (d(i, j)). Suppose the system 
has a unique ground state. Then if the system is gapped, 
so that D(E) = for E < AE, the connected expectation 
value \{Ou\) s — (0) s {u\) s \ is exponentially decaying in I, 
where (...) s denotes the ground state expectation value 
with Hamiltonian TL S . On the other hand, if D(E) cx E a , 
then \(Oul) s — (0) s (ul) s \ is bounded by some constant 
times d{i,0)- a / 2 . 



Definition of Quasi- Adiabatic Continuation 

We introduce the unitary operator 

V(s) (17) 
= S' cxp 1^- ds' dTe-( T/t * )2/2 [ii+(iT) - h.c.]j, 

where the symbol S' denotes that the exponential is S'- 
ordered, in analogy to the usual time ordered or path 
ordered exponentials. We define u s = d s H s = J2i u l> 
and define ufiir) following Ref. [27]: for any operator 
A 

A(t)^A(t) c M-{t/t q ) 2 l% (18) 

A± ^ = hI dtA ^±ir-r (19) 

The time t q will be chosen later. For comparison with 
previous work[27, 28], t q is what was previously called 
yfq/AE; in this work, we also consider the possibil- 
ity of gapless theories where there is no scale AE. 
The time evolution of operators is defined by A(t) = 
cxp[i7i s 't]y4 cxp[— iHs't]. The Hcrmitian conjugate in 
Eq. (17) of ufiir) is uj(— it), and V(s) is a unitary 
operator. 

At t q = 00, the operator V(s) becomes equal to 
V(s) — J2 a *a,s)(*iow,s- To see this, note that at 
t q = 00, we have A + (ir) = A + (ir), where A + (ir) is 
the positive energy part of A taken at imaginary time 
it. That is, in a basis of eigenstates of H s with energies 
E a ,Eb, we have matrix elements A + (ir) a b = A a t,0(E a — 
Eb) exp[— r(E a — £*)], where O(x) is the step function. 
Then, -J™dT[u+(iT)-u-(-iT)} = -{E a - E^dsHs, 
which gives the result of linear perturbation theory for 
the adiabatic evolution of quantum states with a change 
in the Hamiltonian. We instead keep t q finite to define 
a "quasi-adiabatic" evolution, which will transform local 
operators into local operators. To show that keeping t q 
finite maintains the locality we will rely on finite group 
velocity results, while we will use the gap to show that 
we get only a small error in the ground state expectation 
values by taking a finite t q ; detailed proofs of this are in 
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the Appendix, while the physical discussion is given in 
the next section. We will relate the time t q to the scale 
I by Eq. (21) below. 

Results 

For any operator O, we define O(s) = V(s)OV(s)^ , 
where 0(s) has been "smeared out" over a scale I given 
by Eq. (21) below. For a gapped theory we only need 
to take a length I of order the correlation length £ to 
get a small error in Eq. (15). To understand how this 
works, define Q(s) = V(s)V(sy . Then, for a state 
|*iow. s ) in a low energy subspace, (*i O w,s|0(s)|* loWjS ) = 
(*iow,o|Q(s)OQ+(s)|* low , ). If we can show that Q(s) 
has nearly vanishing matrix elements between the low 
energy and high energy states, we can bound the differ- 
ence in expectation values in Eq. (15). The gap between 
the states in the low energy subspace and the high energy 
subspace enables us to show this; for a given gap AE, if 
we perform the quasi-adiabatic continuation sufficiently 
slowly, we can show that (loosely speaking) the matrix 
elements of Q(s) between the low energy and high energy 
states almost vanishes and Q(s) almost acts within the 
low energy subspace as an unitary operator. 

However, there is one complication: in a macroscopic 
system, Q(s) produces excitations throughout the sam- 
ple. Thus, for infinitesimal s, Q(s) = 1 + s Y], e l + 
where the anti-Hcrmitian operator e 1 creates local exci- 
tations near site i. The probability that Q(s) produces 
some excitation acting on the ground state diverges with 
the system size. Strictly speaking, the matrix elements 
of Q(s) between the low energy and high energy states 
are not small. However, the terms e 1 in Q with i suffi- 
ciently far from O can be commuted through O and do 
not affect the final expectation value. We will show that 
the error from the terms e 1 with i near O is exponentially 
small in I divided by the correlation length. The proof 
in the appendix follows this argument, using a triangle 
inequality to make the bound precise. 

If D(E) is not gapped, but has a power law behavior, 
D{E) oc E a for small a, then in a d-dimensional system 



we can still find a local O(s) if a/2 > d + 1 according 
to Eq. (16). This requirement on the exponent can be 
physically understood as follows: if we keep t q finite so 
that O(s) is spread out over a length scale I under the 
quasi-adiabatic evolution, then we have to worry about 
any error in Q(s) on the length scale I. The correlation 
function of u s with O decays as Z~ Q / 2 , and this also is 
how the error terms decay, while the space-time volume 
at scale I is of order l d+1 . Thus, this requirement can be 
understood in terms of the relevance or irrelevance of u s 
at large scales in the given correlation function. 

This physical description is based on two technical re- 
sults. First, we claim that 0(s) is local up to a scale /. 
Specifically, we claim that for any operator Oj which acts 
only on site j, 

W^owwKXtWOjWWow x (20) 

max(exphd(j,0)/ec],cxp[-(d(j,0)/? 9 ) 2 /2]), 

for some constant X\ of order unity and constant l q of 
order t q /c\, where c\ is a characteristic inverse velocity 
of the system. The length £c is a microscopic length of 
order the interaction range of the Hamiltonian and is de- 
fined later. Eq. (20) implies that, while O(s) does involve 
sites more than a distance l q from So, the commutator 
becomes exponentially small once d{j, O) becomes larger 
than l q . To relate l q and /, we ask for the commutator 
to be smaller than the error in Eq. (16) or Eq. (15) once 
d(j,0) becomes larger than I. For a gapless theory, we 
only need to take d(j, O) logarithmically larger than l q 
before the exponential decay becomes much smaller than 
^d+i-a/2. j n j.jjjg casc ^ i j s f orc i er i q _ p or a gapped 

theory, we need that exp[— (l/l q ) 2 } ~ exp[— so that 

l-q/Z-faE/a. (21) 

The correlation length £ is at most of order l/(c\AE). 

If one prefers to have an operator which involves only 
sites within a distance I of So, and exactly commutes with 
sites more than a distance I from So, one may define an 
operator O tr unc(s) such that 



\\Otrunc{s) - 0{s)\\ < X 2 ]T max(cxp [-d(j, O) fa] , exp[- (d(j, 0)/l q ) 2 /2] ) , (22) 

I 



for some constant X2 of order unity. In order for the 
error in Eq. (22) to be of order the error in Eq. (16) or 
Eq. (15), we need to pick I as above: I ~ lq/£ in the 
gapped case, up to logarithmic corrections. 

Second, we claim that O(s) has almost the same ex- 



pectation value as O a diab(s), up to the unitary matrix 
Qo- Specifically, we show that for < s < 1 

|(*low, S |Qo(s) t O( S )g o (s) - Oadiab|*low, S )| (23) 

<2||0||[c 2 ( S ) + C3 ( S )], 
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where c 2 (s),C3(s) are given by Eqs. (32,25) in the Ap- 
pendix. One may verify from the calculation in the 
Appendix that the error term, c 2 (s) + c 3 (s), gives the 
error described above in Eqs (15,16). Note that if 
O = Oi0 2 ...O n , then 0{s) = 1 (s)0 2 (s)...O n (s), so 
the expectation values of products of operators are 
also approximately preserved under this quasi-adiabatic 
evolution. Thus, |(*iow,s|<3o(s) t Oi(s)0 2 (s)...Q (s) - 
,adiab - ■• |*iow, B )| <2||0 1 ||||0 2 ||...[ C2 ( s ) + c 3 ( S )], 
where here the c 2 (s), cs(s) are the error terms appropri- 
ate for the product operator 0\0 2 ... 

EXAMPLES 

We illustrate the quasi-adiabatic continuation by a se- 
ries of examples. We start with local operators, consider- 
ing a system with Z 2 (Ising) symmetry, a fractional quan- 
tum Hall system, and then a system with U(l) symmetry. 
We then repeat the process in the case of non-local, string 
operators. 

Local Ising Model 

The first and simplest example is a quantum Ising fer- 
romagnct in a transverse field. Let the Hamiltonian be 
H = J X)<i j> °f °f + Tlii a i where each site has a spin- 
1/2 and the ferromagnetic interaction J couples nearest 
neighbor spins on the lattice. 

For B = 0, the system has two exact ground states, 
one state which we denote ^ with all spins up and one 
state with all spins down. The gap to the lowest 
excited state above these two states is 2Jq where q is 
the coordination number of the lattice; this state has one 
flipped spin. We define symmetric and anti-symmetric 
combinations ^> S ,A = (1/\/2)(*t ± Tnc states ^ s,A 
are eigenvectors of the operator FJ i erf , with eigenvalues 
±1. This operator, which flips the spin on every site, 
commutes with the Hamiltonian for all B. 

We now consider the quasi-adiabatic continuation with 
K s = -/V , , ofoj + s-B£\of . At s = 0, this has 
the exact ground states ^s,A- For sB small enough, 
the two lowest eigenstates of H s are adiabatic continua- 
tions of ^s,a- Hence, the matrix element of the operator 
d s TL s = B of between these two states vanishes for 
all B, since d s H s commutes with FJ i of and these two 
states are eigenstates of FJi of with different eigenvalues. 

Above these two states, there is a gap to the rest of 
the spectrum. It is known that for small enough B, the 
gap will remain open for all s with < s < 1. Thus, we 
can perform the continuation. The vanishing of the ma- 
trix elements of d s H s between the two low energy states 
implies that Qo(s) is equal to the identity matrix. Wc 
consider the continuation of the operator of . The ground 
state for B = may be taken to be For B = 0, 



we have (\E's|of i of 2 ...o'fJ , I's) = 1 for n even and —1 for 
n odd. The quasi-adiabatic continuation of — > &f(s) 
spreads out the of operators over a distance I. Cer- 
tainly I is less than the linear size of the system. Then 
(*s, s |of i (s)of 2 (s)...of n (s)|<I's iiS ), again equal to 1 or -1, 
up to some error, depending on whether n is even or odd. 
Here, = V(s)^s- The error is exponentially small 
in l/£. We conjecture that the operators of (s) corre- 
spond to block spin operators: they are equal to plus or 
minus 1 depending on whether the average spin over a 
correlation volume is positive or negative. 

If we instead consider the Hamiltonian H s = 
JJ2<i, 3 > a i a j + sB J2i u h tne operator d s H s will have 
non- vanishing matrix elements between the states ^>s,A, 
and thus we must start with states ^fj^j. to per- 
form the continuation if we want to have Qq(s) be 
equal to the identity. In this case, of (s) = of and 
(* T;S |o-?(s)cr| 2 (s)...o| ! n (s)|* T;S ) = (+1)™, where we have 
assumed that the sign of B is such that ^>^ s is the ground 
state, rather than ^i, s - For the opposite sign of B, the 
correlation function is instead (—1)". 

The two cases, depending on the different ways to add 
the magnetic field transverse or parallel to the z-axis, 
lead to different ground state correlation functions. In 
the second case, the correlation function of the continued 
operator is not a very interesting result; the unit opera- 
tor would have the same correlation function. However, 
in the first case, the ability to find a continuation of the 
operator with the given correlation functions is a much 
more interesting result. There are long-range correla- 
tions in the operator, since (\I/s >s |of (s)of (s)|5 , 5. s ) — 
(*s,<,K(*)|*sr,«X*fir,.|<T?(a)|*sf >a ) is non-vanishing 
even when sites i\,i2 are far from each other, and in 
particular even when the distance between ii,i 2 much 
larger than I. This implies [27] the presence of another 
state close in energy to the ground state when the 
magnetic field is added in the transverse direction. 

However, we would like to show the ground state de- 
generacy in another way, based directly on continuing 
the states ^ s.a- This way will be very valuable for more 
complicated systems such as the quantum Hall system. 
We will consider a more general class of Hamiltonians 
TL S : we consider arbitrary Hamiltonians which are sums 
of local terms, and which commute with TJ i of. We show 
that, under the assumption that the gap between the 
two lowest states and the rest of the spectrum remains 
open, the energy difference between the continuation of 
the two lowest states is exponentially small. As long 
as that gap remains open, it is possible to adiabatically 
continue these two states, giving states ^s. s = y(s)*s 
and ^a,s = V(s)^a as eigenstates of 7i s - Here, we rely 
on the fact that Ti. s commutes with Yli &\ and thus has 
vanishing matrix elements between the states. 

To compute the difference in energies for given 
s, we compute (*a, s |W s |^, s ) - (*5. S |W S |* S . S ) = 
(y A \V(s)iH s V(s)\* A ) - (*s\V{s)iH s V(s)\* s )- Thus, 
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V(s)^H s V(s) defines the continuation of the operator 
Ti s from s back to and V(syH s V(s) defines a quasi- 
adiabatic continuation from s back to 0: (^s.sl^sl^s.s) 
is approximately equal to s\V (s)^H s V (s)\$> s ) ■ The 
error in this continuation is exponentially small in Z/£. 
The operator H s is a sum of local operators, while the 
operator V(syH s V(s) is a sum of terms spread out 
over length scale I. The only operators O such that 
(*a|0|*a> - (*s|0|*s> ^ are operators which flip 
every spin in the system and thus have non-vanishing 
matrix elements between and "J^. Thus all lo- 
cal operators have vanishing matrix elements between 
the two states "I^ and "J^. In particular, if the 
length scale I is smaller than the system size L, then 
(*A\V(8)iH a V(a)\* A ) = (*s\V(8)iH a V(8)\*s)- Thus 
we can pick I just smaller then the system size to show 

that (*A\V(8)iH B V{8)\* A ) - (*s\V{8)lH a V(8)\* s ) is 
of order ||W S || exp(— L/£) ~ L d cxp(— £/£). Here, the 
bound for this system is not a very tight bound: one 
expects the level splitting to be exponentially small in 
(L/C) d instead. 

The key steps in this argument were that (1): matrix 
elements of operators which commute with FJ i of vanish 
between and ^a] and (2): all local operators have 
the same expectation values in state ^a as in state *s. 

Quantum Hall Effect 

We now turn to the case of the fractional quantum 
Hall effect. Consider a system with no disorder on a 
torus at filling factor p/q, with p and q coprime. The 
magnetic translation group implies at least a g-fold de- 
generacy of the ground state[18]. Assume that in the 
absence of disorder there is a g-fold degenerate ground 
state, with a gap to the rest of the spectrum. Now, 
consider adding disorder to the system, defining Ti s = 
Ho + s J dxdyU(x,y)^(x,y)^(x,y), where U(x,y) is a 
disorder potential and Ho is the clean Hamiltonian. Wen 
and Niu argued[18] that to first order in sll(x,y) the 
splitting between the q-fold degenerate states was expo- 
nentially small. We will use the continuation to show that 
the splitting is exponentially small for non- vanishing dis- 
order strength under the sole assumptions that the gap to 
the rest of the spectrum remains open and that at s = 
all local operators have the same expectation value up 
to exponentially small terms in the q lowest states and 
that at s — all local operators have exponentially small 
matrix elements between the q lowest states. We note 
that since these last statements involve only s — they 
can be checked in specific model systems without disor- 
der, and in fact were checked in [18] when they showed 
the exponentially small splitting of the q lowest states in 
linear perturbation theory. 

Suppose, then that for < s < 1 there are q 
states, for n = 0...q — 1, and a gap to the rest of 



the spectrum. We wish to show that these states are 
close in energy. Thus, we compute (^ m ,s\1~ts\^m,s) — 
(^n^lHsl^n^), for some m,n = 0...q — 1. As above, 
V(s)'' i HsV(s) defines the continuation of the operator 
H s from s back to 0, and V(s)^H s V(s) defines a quasi- 
adiabatic continuation from s back to 0. 

However, unlike the case of the quantum Ising model, 
we do not have any symmetries to make matrix el- 
ements of H s and d s H s vanish between the low ly- 
ing states. Thus, we do not have any control on 
the matrix Qo(s). Then, (^m.slHsl^m^) is equal to 
(^ m ,o\Qo(sYV(syH s V (s)Q (s)|* m ,o), up to an error of 
order | \H S \ | exp(— lj £) . Since we are continuing from non- 
zero a to s = 0, now the matrix Qo acts within the low 
energy sector of Ho ■ 

However, this matrix Qo(s) causes no problem. As 
before, for I < L, V(s)^H s V(s) is a local operator, and 
then under the assumptions above, the expectation value 
(^ m ,o\Qo(sYV(sYHsV( S )Qo(s)\^ m ,o) is independent of 
Qo(s), up to exponentially small corrections, thus giving 
the desired result. 



Local Rotor Model 

The next example is a system with a continuous sym- 
metry. We take the Hamiltonian H s = k(s)~ 1 J2i^i + 
k(s)^2 <i zizj, where z% is a continuous complex field 
with \zi\ — 1, and Ilj is a momentum with [zi,IL;] = iz%. 
The parameter k(s) is an s-dependent stiffness of the field 
z. 

We pick a quasi-adiabatic continuation using k(s) = 
ko(ki/ko) s , so k(0) = ko and fc(l) = k\, thus d s H s = 
Hh/kon-Hs)- 1 n? + k(a) E< lJ > zizj]. We choose 
the initial ko to be large compared to the final k\. We 
assume that we k\ is sufficiently large that the system 
is still in a phase with gapless excitations and algebraic 
correlations. In this phase, we can compute the density 
of states by writing z = exp(i^) for some 4> with a Gaus- 
sian action for <f>. At low energy, the system acquires 
rclativistic invariance, and thus in d dimensions the den- 
sity of single particle states at energy E is of order E d ~ l . 
The matrix element of <f> between the ground state and 
such a single particle state is of order E^ 1 ' 2 , and thus the 
integrated density of states created by <j) is D{E) oc E d ~ l . 
The integrated density of states below energy E created 
by dsHs instead is KD(E) - ln(fc / h) E 2d + 2 . Then, 
a/2 = d + 1 and we are in a marginal case: the error in 
the continuation is of order log(L/Z). 

Emergent Discrete Gauge Theories 

We now turn from these theories with local operators 
to emergent gauge theories. We consider the Hamiltonian 
of the emergent Zi gauge theory in the introduction. 



11 



For Ji = J2 = 0, we have 4 exactly degenerate ground 
states on a torus, and a gap to the rest of the spectrum. 
For non-zero Ji, J2, the deformed Hamiltonian is still lo- 
cal 7i s = X^i^s) where HI is a local operator defined 
near the site i. We can use the same reasoning used in 
FQH states to show that Tig still has 4 exactly degener- 
ate ground states on a torus. However, here we will use 
a slightly different approach. Continuing 7i s from s back 
to 0, we find that 

(* m , s |W:i*„, s ) = (* m ,o|V t (a)W*V(*)|tt n , >, 

where \^ n ,s), n — 1, 4, are the 4 low lying states of 
H s . Due to the energy gap separating the 4 low lying 
states with rest of states, the operator V' f (s)HlV(s) is 
almost a local operator. More precisely, up to an error of 
order e~ 1 ^ and a unitary rotation Qo(s) between the 4 
low lying states \® n ,o), VHs)H s V{s) can be replaced by 
a truncated operator H % which only acts on sites within 
a distance I from the site i. Thus we have 

(*m,s\K\*n,s) = (* m ,o|Oo(s) t J tf i OoW|*„,0>+0(e-^). 

Since Qo(s)\^ n ,o) are the ground states of the exactly 
soluble model (1), they from a irreducible representa- 
tion of the algebra of the large-closed-string operators 
(8). We can choose the length I over which the opera- 
tors are smeared to be one quarter of the linear size L of 
the system. In this case H % will be local enough that we 
can choose the positions large-closed-string operators to 
avoid any overlap between H % and the closed-string oper- 
ators. So ft 1 commutes with the closed-string operators. 
As a result, H l must be proportional to the identity oper- 
ator within the irreducible representation. This way, we 
have shown that (^ r miS |Wg|'J' njS ) oc S mn up to an error 
e -V 4 £ ; which implies that (\&m, s |W s |^ n , s ) ^ &mn U P to 
an error L 2 e~ L / 4 ^. So the energy splitting between the 
4 low lying states is less than L 2 e~ L / 4 ^ for the deformed 
Hamiltonian 7i s . 

We can also continue any closed-string or dual closed- 
string operators iS(C c i OS ed) and obtain 

(*m,o|-S(Cclosed)| *n,o) = (*m, s \V(s)S(C closed )Vl(s)\H> niS ), 

Then there exist dressed closed string operators 
S' dre (C c io S od) that have a width I such that 

(* m , s |go(s) t 5 dro (c closod )g ( s )|*„, s > 

= <* m ,a|^(*)S(C c l OB ed)^ t («)|*n,a> + 0(e~ 1 ^), 

Here Qa(s) is a unitary rotation between the 4 low ly- 
ing states \^ n ,s) and Qo(s) is independent of the closed- 
string operators. This implies that the dressed string and 
dual string operators £ dre (C c iosed) have the same algebra 
among the low energy states of the perturbed Hamilto- 
nian as the original operators S'(Cciosed) do for the orig- 
inal Hamiltonian (up to an error exp(— //£))• 



The above also implies that the expectation of the 
dressed closed-string operators in the ground state of the 
perturbed Hamiltonian satisfies the zero law (3), up to 
an error of order \S\ exp(— where \S\ is the string 
length. We see that the error is exponentially small for 
long strings since / can be chosen to be a fraction of \S\. 

Our final task is to show the invariance of certain states 
under a deformed gauge transformation for J 2 7^ 0. First 
consider J\, J2 = 0. The Hamiltonian has the 4 ground 
states which are invariant under the local gauge trans- 
formation Wi given by Eq. (2). We can create gauge in- 
variant excited states by acting on the 4 ground states by 
open dual string operators. Following Eq. (5), we define 
an open dual string operator as S(C) — J[ i cross g of , 
where now the dual string C is open, with two end- 
points. The states created by acting with S(C) are still 
invariant under the local Z2 gauge transformations, since 
[S(C), Wi] = 0. They are excited states which introduce 
Z2 gauge flux at the endpoints of the open string. By 
acting on the ground state with the open dual strings, 
we can create all gauge invariant states, and for U 3> <?, 
these are the lowest energy excited states. Continuing to 
non-zero J\ at J 2 = still leaves these states invariant 
under the local Z 2 gauge transformation. Continuing to 
non-zero J 2 breaks the local Z2 gauge invariance. How- 
ever, these states are invariant under the deformed local 
Zi gauge transformation Wj{s). To see this, we use the 
continuation to show that (^iow,s|l / I / i"(s)| x I , iow,s) = 1, up 
to exponentially small error in Since Wj(s) is uni- 
tary, this implies the gauge invariance of the ground state 
under the deformed gauge transformation up to exponen- 
tially small error. The gauge invariance, up to the same 
small error, of the states created by acting on the ground 
state with operators S(C,s) = V (s) S (C , s)V (s)t then 
follows from the exact commutator [S(C, s), Wi(s)] = 0. 
The size of the deformed gauge group is, however, much 
smaller than the original: if the linear system size is 
L, there are (L/l) d different local gauge transformations 
rather than L d as is the case at J2 = 0. 

Emergent Continuous Gauge Theories 

Our final problem is the theory with the U(l) gauge 
symmetry, Eq (13). The absence of a gap makes it much 
more difficult to obtain results on this system. We will 
obtain only one result, the existence of a deformed U(l) 
gauge invariance of the system. 

We consider some gauge transformation acting on a 
site, Wi^ = e l ^ Ql . When J 2 = 0, the ground state 
has an exact gauge invariance: (\&o|Wi,4>|\&o) = 1. If 
we continue to non-zero J2, the ground state breaks the 
gauge invariance. However, we claim that the expecta- 
tion value of the continued operator, ( x I'o,s|W^/,0(s)|*o,s)) 
is still close to unity. Unlike the Z 2 case, this does not fol- 
low simply from the results derived previously as there 
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is no gap, and the low energy density of states is too 
high to use the results for gapless systems. However, 
it is still possible to show this result. We only very 
briefly sketch the argument, leaving a more detailed pre- 
sentation for future work. The idea is as follows: at 
J2 = 0, the system has gauge non-invariant states at an 
energy of order U above the ground state. Under the 
continuation, at finite t q , the operator V(s) will take the 
ground state at s = into some state which is a mix of 
ground and excited states. That is, the expectation value 
(*o,.|W/,*(«)|*o,a> = (*o[Q(s) t ^ / , Q( S )|* o ), where 
the unitary matrix Q(s) — V(s)W(s) mixes the ground 
state with the excited states. However, since all the low- 
lying excited states at J 2 = are gauge-invariant, it is 
no longer necessary to show that Q(s)^Wi i< j,Q(s) is close 
to Wi^\ it suffices to show that we do not mix in the 
states which are not gauge invariant under Wj^ and all 
such states are at energy of order U above the ground 
state. At infinitesimal s, the mixing into such states is 
exponentially small, since the energy U acts like a gap; 
however, one has some mixing into the low-lying gauge- 
invariant states. As s increases, one mixes into progres- 
sively higher energy states, until eventually one begins 
to excite the gauge non-invariant states. However, if t q 
is sufficiently big (roughly of order U -1 ) the mixing into 
the gauge non-invariant states can be bounded at s = 1, 
thus obtaining the desired result. 

Now consider the excited states. At J2 = 0, the gauge 
invariant excited states of the system can be obtained 
by acting on the ground state with operators of the 
form O — Yli e lc ^ iLi , where the phase fa is some arbi- 
trary function of the leg. These operators commute with 
Wj^ and therefore at non-zero J2 the continued opera- 
tors commute with the continued gauge transformation: 
[0(s),Wi v< p(s)] = 0. Therefore, there are a class of ex- 
cited states for non-zero J2, namely the excited states 
created by acting on the ground state with the operators 
0(s), which are also gauge invariant under the deformed 
gauge transformation, up to the same error. 

As in the Z2 case, the size of the deformed gauge group 
is much smaller than the original gauge group. The orig- 
inal gauge group had L d different generators, while the 
deformed group has only (L/l) d such generators. A more 
careful analysis should be able to then use this deformed 
gauge invariance to show that the gapless photon is pro- 
tected. This is also a job for the future. 

The compactness of the gauge group was important 
here. Consider a non-compact U(l) theory, with La- 
grangian(l/2)(a AJ A,-^A Al ) 2 + (A/2)(a AJ ^)2+MA AJ ^. 
The term in A is a gauge fixing term and M is a term 
which breaks the gauge symmetry. With this quadratic 
action, it is easy to verify that a non-zero M opens a gap. 
Why doesn't the continuation work here? The reason is 
that this theory has no gap to the gauge non-invariant 
states, unlike the compact cases before, despite the gauge 
fixing term. The gauge fixing term does not open a gap to 



the gauge non-invariant states in this case, instead it adds 
gapless longitudinal and scalar photons to the theory. To 
say it differently, in a compact theory, the charge is quan- 
tized and the open string states have a minimum possible 
energy because they necessarily terminate in an end with 
a charge that is a multiple of the charge quantum. In a 
non-compact theory, the charge may be arbitrarily small. 
So the energy cost is also arbitrarily small. 



DISCUSSION 



The standard wisdom is that as long as gaps remain 
open, a system does not have a quantum phase transi- 
tion and thus the long-distance structure of correlation 
functions remains the same. We have shown a precise 
form of this statement. We have found that, by appro- 
priately dressing operators, the long-distance structure 
can in fact be preserved to a much greater degree than 
one might have expected. In particular, we have shown 
the presence of the zero law (3) for gauge theories in 
the deconfined phase, and we propose this as a test of 
confinement. Further, we have considered the stability 
of topological order under perturbations of the Hamilto- 
nian, and shown that the order is robust unless the gap 
to the rest of the spectrum (local excitations) closes. 

Topologically ordered states are described by emergent 
gauge theories at low energies. The topological order is 
closely tied to the emergent gauge invariance of the low 
energy gauge theories. From the point of view of the low 
energy gauge theory, our result shows that the emergent 
low energy gauge invariance is topological. It cannot be 
broken by any local perturbations in the parent bosonic 
model. We hope this result will shed light on the true 
meaning of gauge invariance and gauge theory. 

The continuation is also useful for systems without 
emergent local gauge structure, and may have a wider 
applicability. For example, an outstanding question 
is to prove that in some neighborhood of the AKLT 
point[29, 30], a spin-1 chain remains gapped[31]. The 
continuation might be useful in doing this, or at least 
in showing, under the assumption of the existence of a 
gap, the persistence of string order [32] throughout the 
Haldane phase. 
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PROOF OF LOCALITY AND APPROXIMATION 
RESULTS 

Locality Result 

To show Eq. (20), wc use the finite group veloc- 
ity result, proven in [26, 28]. This result uses the 
finite-range conditions on the Hamiltonian above to 
bound the commutator ||[^4(i), 5(0)]||, where A(t) = 
cxp(iH s t)Aexp(—iH s t). One can show that this com- 
mutator is exponentially small for times t less than c\l 
where / is the distance between A and B and c\ is some 
characteristic inverse velocity which depends on J, R, 
and the lattice structure. The specific bound is that 
1 1 [A(t), B(0)\ | < 1 1 A\ 1 1 |B| | Y,j9{t, d{A, j)), where the sum 
ranges over sites j which appear in operator B and where 
the function g has the property that for \t\ < c\l, g{c\l, I) 
is exponentially decaying in / for large / with decay length 
£c for some constant £c which is of order R. Recall that 
d(A,j) is the minimum over sites i acted on by A of the 
distance d(i, j). 

Before giving the proof of Eq. (20), we give a phys- 
ical description. The finite group velocity result has a 
very simple interpretation. Consider a local operator A. 
Under time evolution, we get an operator A(t) which 
"spreads out" over space as time passes. The finite group 
velocity result implies that for finite t the operator A is 
still local up to some length t/c\ (here c^ 1 is some char- 
acteristic velocity of the system) in the following sense: 
the commutator of A(t) with any operator B is expo- 
nentially small if B is at least distance t/c\ from A(0). 
This applies in particular for A — u s . Thus, in the def- 
inition (17) of U, the operators wj(ir) are local in the 
same sense: Eq. (19) gives wJ(«t) as an integral over t 
of u s i (t) and for t » t q the integral is cut off exponen- 
tially, while for t ~ t q the u a > (t) are local up to length of 
order t q /c\. Then, we define O(s) by the unitary trans- 
formation V(s). We view this unitary transformation as 
defining a fictitious time evolution with time parameter 
s and Hamiltonian given by the exponent of Eq. (17). 
We have just established that this exponent is local and 
we can then apply the finite group velocity result to this 
evolution to show that O(s) is also local. The rest of this 
subsection consists of a few precise error bounds following 
these statements. 

We compute the commutator [m* + (it), Oj] where Oj 
is some operator which acts only on site j, and 
where w* + (ir) is defined following Eq. (19) taking 
A = u\. We separate the integral over times t 
into times with t < c\l, where I = d(i,j) — R, 
and times with \t\ > c\l. For the first set of 
times, we have (2n)^ 1 f, t , <Cil dt(it + T)\\[u s (t), Oj}\\ < 
(27r)- 1 ||u 4 s ||||0||exp(-?/^ c ). For the second set of 
times, we have (2tt)^ 1 j dt(it + T)\\[u s (t),Oj]\\ < 
K~ 1 \K\\\\0\\(V2^t q /c 1 l)exp[~(c 1 l/t q ) 2 /2]. Thus, for 



large /, we find ||[w* + (it), Oj]\ \ is exponentially decaying 
in d(i,j) with decay length £c- Note that exp(— l/£c) ~ 
cxp[-( Cl Z/g72] for I « 2(t 9 / Cl )7£c- 

Now, here is the trick. We regard Eq. (17) as 
defining the "time" evolution of states, where the pa- 
rameter s' is an effective time parameter and D = 
i Jq 00 drexp[— (t /t q ) 2 /2][d^,(iT) — h.c] is some effec- 
tive s'-dependent "Hamiltonian", so that V(s) = 
<S'exp[— i f° ds'D]. We write D — J2i^\ where 
D l = / ° c drexph(T/t 9 ) 2 /2][u 4 s +(iT) - h.c.]. Then, 
||[£> 4 ,0 3 ]|| < \\Oj\\F(d(i,j)), where the function F(l) is 
exponentially decaying as exp[— l/£c] for large I and de- 
caying as (\Z2nt q /cil) exp[— (al/t q ) 2 /2] for small I. 

This exponential decay is in fact good enough to prove 
the finite group velocity result [26] using D as an effec- 
tive Hamiltonian; it is not necessary that D l act only on 
sites within some finite range, but an exponential decay 
also suffices. Following [28], define Gi by the differential 
equations, for s > 0, d s Gi(s) = J2j Gj(s)F(d(i,j)) with 
initial conditions Gj(0) = 2 if i e So and G,(0) = 
otherwise. Then, one can show that ||[0(s), Oj]\\ < 
\\0\\\\Oj\\Gj(s). Solving the equations for Gj, one ar- 
rives at the bound (20). 

One can define the operator O tr unc by setting 
Vtrunc(s) = 5'cxp[-z/ s ds'Tr : ,^( Oj )> ( (L))], where the 
trace is the trace of operator D over sites j with d(0, j) > 
I. Then, set O trU nc{s) = Vtrunc{s)OV t runc{s) f , getting 
Eq. (22). 

Approximation Result 

The idea behind the proof of the approximation re- 
sult is that, for any site i, the difference between u 1 +{it) 
and ul + (ir) can be made small by taking large enough 
t q . Here, u l s + is the positive energy part of u s and w* + 
is defined following Eq. (19) with A — u\. The differ- 
ence between these two is closely related to e\ as given 
below. Eq. (17) involves summing over all sites i, but 
sites i which are sufficiently far from O will turn out to 
have little effect on defining O(s). Thus, the task in this 
subsection is to figure out how much difference there is 
between u l s + (ir) and w* + (ir), and then sum that error 
over sites near enough to O, giving the difference be- 
tween the quasi-adiabatic continuation and the adiabatic 
continuation. Since the adiabatic continuation preserves 
expectation values, this will give an estimate in the error 
in the expectation values. We now do this carefully 

The proof of Eq. (23) involves defining an additional 
operator and using triangle inequalities. We define 

Vtis) = S'exp{- / ds' / dr^cxp[-(r/t 9 J 2 /2] ><24) 
Jo Jo j 

[uit(iT)-h.C.]}, 

where now t qi may depend on i and we define u l s + {it) by 
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u\(t) = u\ (i) exp[— (t/t qi ) 2 /2], again using the t Qi which 
depend on i. 

We then pick t Qi = t q for d(i, O) < 21, and t Qi = t g + 
ci[d(i,0) - 21} otherwise. Define Oj(s) = Vi(s)OVi(s)^ . 
Thus, V/(s) has a t qi which increases the further one gets 
from operator O. While the operator V(s) would create 
local excitations everywhere acting on a ground state, 
Vi(s) only creates local excitations near O. 

Using a triangle inequality, | (*i w,s|Qo(s)0(s)Q (s) t - 

O ad iab(s)|*low, S )| < |(*low, S |Oo(s)0(s)Q (s) t - 

Qo(s)0 ; (s)Q (s) t |*iow,s)| + |<*iow,.|Qo(*)Oi(*)Qo(«) t - 
Oadiab(s)|*iow,s)|- Tne difference between V t (s) 
and V(s) is the excitations far from O. That is, 
if for small s U(s)U(s) t = 1 + sJ2 i e l + ... and 
Vi(s)V(sy = l + s^ej + then ej and e l differ 
only for i far from O. However, for i far from O, 
the e 1 commute through O (as discussed physically 
before) and so it is possible to bound the difference 

K*low,B|Qo(«)0(*)Qo(*) t - Oo(«)Ol(*)Oo(*) t |*low,B>|. 

Precisely, to bound the first difference, we note that the 
difference between the definition of Oi(s) and 0(s) has 
to do terms u l 8 with sites i which are at least a distance 
21 from O. Using the locality bound, one can bound the 
commutator of O(s) and Oi(s) with u\ for these sites. 
This gives \ \0(s) — Oi(s)\ \ < c 3 (s), where 

c 3 (s)=X 3 max(exp[-d0',O)/fc], (25) 

j,d{i,j)>2l 

exp[-(d(i,0)// g .) 2 /2]), 

for some constant X 3 and where l gj is of order C\t qj . 
We now bound the difference 

|(*low,s|Qo(s)0(( s )Qo(s) t - Oadiab(s)|*low, S )|- To do 

this, it suffices to bound |(*iow,s<9o(s)V;(s)-(*iow,s^(s)|- 
This is equal to | (*iow L oU(s)tQ ( s )VKs) - 
(*iow,o|- The operator Vj(s) is equal to 

S'exp{- f°ds'(d s ,Vi(s'))Vi(s'y}. This equals 

V{s)S t exp{-fids'V{8 t MdsV l {8 t ))V l {8'y 
(d s ,V(s'))V(sy]V(s)}. Thus, 

V(s)iQo(8)Vt(s) (26) 
= [V(s)iQ (s)V(s)] x 

S' exp{- J' d8'V{ayei( 8 ')V(a')}, 

where 

ei (s') = (d a ,V l (s'))V l (s')W(s') - (d s ,V(s'))V(s')l (27) 

We have grouped the operators V(s)^Qo(s)V(s) to- 
gether in (26) for a reason: the matrix Qq(s) is an op- 
erator between the low energy states of H s so therefore 
the operator V(s)^Qo(s)V(s) is an operator between the 
low energy states of Wo- Now we turn to the exponen- 
tial S'exp{- f s ds'V(s'yei(s')V(s')}. We want to show 



that this operator is also equal to, up to some bounded 
error, an operator between the low energy states of Ho- 
Define Phigh to project onto the high energy states of 
Ho- Then, we can pick Qo(s) such that 

|(*iow,o^(s) t g (s)Vi(s) - (*iow,o| (28) 

< / ds'\(* lowfi V( S yei(s')V( S ')P high \, 
Jo 

We will bound the integral of Eq. (28) below; combining 
this bound with Eq. (25) will give Eq. (23). 
Using linear perturbation theory, 

poo 

(d s/ V(s'))V( S y = - / dr[w+(ir) - h.c] + P, (29) 
Jo 

where P only has non- vanishing matrix elements between 
states of the same energy: P a b = if E a ^ E^. 
The error ej(s) = — P + Y2i e K s ) wri ere 

pOO 

efts) = - / dT[exp[-(T/t qt ) 2 /2]ul+(tr) - «'+(ir) - ft-cpO) 
Jo 

For s = 0, e\{s) = defined above. 

We now bound the projection into the high en- 
ergy sector |(*iow,s'eJ(s')i\i g h|. We can show by 
performing some elementary integrations [27] that, 
for any eigenstate ^/ a ,s with energy E a > i~/t q ., 
|(* low , s |e X p[-(r/^) 2 /2]^ + (iT) - uj+(ir)|* B , B )| < 
cxp[-(rA 9i ) 2 /2]cxp[-(^i? a ) 2 /2]|«)o a |, where \(4) 0a \ 
is the absolute value of the matrix element of u l s be- 
tween state ^low.s and * a ,s- Similarly, for E a < r/t 2 ., 
|(* low , s |exph(T/^) 2 /2]^+(ir) - uj+(iT)|*„, 8 )| < 
cxp[— ri? a ]|(u*)oa|- Integrating over r and summing over 
states outside the sector of ground states, using the 
bound on density of states, we have 

K*iaw, B eks)iW»| 2 < / dEp(E) K31) 

{(K/E)2cxp[-(t qi E) 2 } + Kt qi cxp[-(t qt E) 2 /2]V2^)} 2 . 

Summing over sites i and using Eq. (28) we obtain the 
bound |(*iow,oU(s) t Oo(s)U;(s) - (*i ow ,o| < c 2 (s) where 
we define 

c 2 (s) = S E(/ dEp(E){(K/E)2cM-(t qz E) 2 } + (32) 

Kt qi exp[-(^S) 2 /2]V2^} 2 ) . 
This completes the calculation. 
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